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ABSTRACT 

We present an absolute parallax and relative proper motion for the fundamental distance 
scale calibrator, 8 Cep. We obtain these with astrometric data from FGS 3, a white-light 
interferometer on HST. Utilizing spectrophotometric estimates of the absolute parallaxes of our 
astrometric reference stars and constraining 8 Cep and reference star HD 213307 to belong to 
the same association (Cep OB6, de Zeeuw et al. 1999), we find 7r a (, s = 3.66 ± 0.15 mas. The 
larger than typical astrometric residuals for the nearby astrometric reference star HD 213307 are 
found to satisfy Keplerian motion with P = 1.07 ± 0.02 years, a perturbation and period that 
could be due to a F0V companion ~ 7 mas distant from and ~4 magnitudes fainter than the 
primary. Spectral classifications and VRIJHKT 2 M and DD051 photometry of the astrometric 
reference frame surrounding 8 Cep indicate that field extinction is high and variable along this 
line of sight. However the extinction suffered by the reference star nearest (in angular separation 
and distance) to 8 Cep, HD 213307, is lower and nearly the same as for 8 Cep. Correcting 
for color differences, we find <Ay > = 0.23 ± 0.03 for <5 Cep, hence, an absolute magnitude 
My = -3.47 ±0.10. Adopting an average V magnitude, <V> = 15.03 ± 0.03, for Cepheids 
with log P = 0.73 in the LMC from Udalski et al. (1999), we find a V-band distance modulus 
for the LMC, m-M = 18.50 ± 0.13 or, 18.58 ± 0.15, where the latter value results from a highly 
uncertain metallicity correction (Freedman et al. 2001). These agree with our previous RR Lyr 
HST parallax-based determination of the distance modulus of the LMC. 

Subject headings: astrometry — interferometry — stars: distances — stars: individual (8 Cep) 
— stars: binary — distance scale calibration — LMC 
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Introduction 



Many of the methods used to determine the distances to remote galaxies and ultimately the size, 
age, and shape of the Universe itself depend on our knowledge of the distances to local objects. The 
most important of these are the Cepheid variable stars. Considerable effort has gone into determining the 
absolute magnitudes, My, of these objects (see the comprehensive review by Feast, 1999). Given that the 
distances of all local Cepheids, except Polaris, are in excess of 250pc, most of these My determinations 



rely on large- number statistics, for example Groenewegen fc Oudmaijer (2000) . Lanoix, Paturel, fc Gamier 



(1999), and Feast (1997). Gieren et al. (1993) used Cepheid surface brightness to estimate distances and 



absolute magnitudes. For Cepheid variables, these determinations are complicated by dependence on color 
and metallicity. Only recently have relatively high-precision trigonometric parallaxes been available for a 



very few Cepheids (the prototype, S Cep and Polaris) from HIPPARCOS (5 Cep = HIP 110991, [Perryman 
et al. (1997)[ ). We have determined the parallax of S Cep with FGS 3 on Hubble Space Telescope with 
significantly higher precision. Additionally, our extensive investigation of the astrometric reference stars 
provides an independent estimation of the line of sight extinction to S Cep, a significant contributor to the 
uncertainty in its absolute magnitude, My. 

In this paper we briefly discuss (Section ^) data acquisition, analysis, and an improved FGS 3 
calibration; present the results of spectrophotometry of the astrometric reference stars required to correct 
our relative parallax to absolute (Section |J); and derive an absolute parallax for S Cep (Section iA). Finally, 
we calculate an absolute magnitude for <5 Cep (Section |5.3| ) and apply it to derive a distance modulus for 



the Large Magellanic Cloud (Section 5.4) 



Bradley et al. (1991) and Nelan et al. 2001 provide an overview of the FGS 3 instrument and Benedict 



et al. (1999) describe the fringe tracking (POS) mode astrometric capabilities of FGS 3, along with data 
acquisition and reduction strategies used in the present study. We time-tag our data with a modified Julian 
Date, mJD = JD - 2444000.5. 



2. Observations and Data Reduction 

Figure [l] shows the distribution in RA and Dec of the five reference stars and S Cep. Seven sets of data 
were acquired, spanning 2.44 years, for a total of 127 measurements of S Cep and reference stars. Each data 
set required approximately 40 minutes of spacecraft time. The data were reduced and calibrated as detailed 
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Benedict et al. (1999) and McArthur et al.(2001). In a recent paper (Benedict ct al. 2002) we described 



a technique used for these data in which we employ a neutral density filter to relate astrometry of very 
bright targets to faint reference stars. At each epoch we measured reference stars and the target, 5 Cep, 
multiple times. We do this to correct for intra-orbit drift of the type seen in the cross filter calibration data 



in our recent paper reporting a parallax for RR Lyr ( Benedict et al. 2002 , Fig. 1). Data sets 2 and 5 were 
each also afflicted with an episode of non-monotonic drift, possibly due to mechanism (filter wheel) motion. 
Fortunately, the strategy of multiple repeats of the observation sequence within each data set permitted 
generally satisfactory correction. 

All of these observations were acquired with FGS 3, and all before 1998. Since 1999, the prime 
astrometer aboard HST has been FGS lr, installed during the 1999 servicing mission. To calibrate the 
Optical Field Angle Distortion (OFAD) for FGS lr we have, during the last three years, secured over 30 
observation sets of M35, our calibration target field. These new observations have extended our time base to 
over ten years and allowed us to refine the proper motions of calibration stars in M35. This has resulted in a 
more precise star catalog, which in turn has allowed us to improve our FGS 3 OFAD. Applying this revised 



calibration to the Barnard's Star data presented in Benedict et al. (1999) , we find a parallax difference of 
0.07 mas and a proper motion difference of 0.7 mas yr _1 , each a change of about one part in 10,000. We 
use this new calibration for these 5 Cep data. 



3. Spectrophotometry Absolute Parallaxes of the Astrometric Reference Stars 

Because the parallax determined for 8 Cep will be measured with respect to reference frame stars which 
have their own parallaxes, we must either apply a statistically derived correction from relative to absolute 
parallax (van Altena, Lee & Hoffleit 1995, hereafter YPC95) or, preferably, estimate the absolute parallaxes 
of the reference frame stars seen in Figure [j]. With colors, spectral type, and luminosity class for a star one 
can estimate the absolute magnitude, My, and V-band absorption, Ay. The absolute parallax is then, 

n abs = 10 -(v-m v+5 -a v )/5 (1) 

The luminosity class is generally more difficult to determine than the spectral type (temperature class) . 
However, the derived absolute magnitudes are critically dependent on the luminosity class. To confirm 
the luminosity classes we employ the technique used by Majewski et al. (2000) to discriminate between 
giants and dwarfs for stars later than ~ G5, an approach whose theoretical underpinnings are discussed by 
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Paltoglou fc Bell (1994). 



3.1. Photometry 

Our band-passes for reference star photometry include: BVRI, JHK (from preliminary 2MASS[] data), 
and Washington/DDO filters M, 51, and T2 (obtained at McDonald Observatory with the 0.8m Prime 
Focus Camera). We transform the 2MASS JHK to the Bessell (1988) system using the transformations 
provided in Carpenter (2001) . In Tables |l| and || we list the visible, infrared, and Washington/DDO 



photometry for the 5 Cep reference stars, DC-2 through DC-7. DC-2 was too bright for 2MASS and our 
Washington/DDO and RI photometric techniques. 



3.2. Spectroscopy 

The spectra from which we estimated spectral type and luminosity class come from the New Mexico 
State University Apache Point Observatory^. For all but reference star DC-2 classifications were obtained 
by a combination of template matching and line ratios. Differing classifications for reference star DC-2 (HD 



213307) have been reported in the literature. We consider the classifications of both Lutz & Lutz (1977) and 
Savage et al. (1985). Table [| contains a list of the spectral types and luminosity classes for our reference 
stars, rank ordered by estimated distance. We discuss our estimation of the <Ay > in the next subsection. 

In Figure |^ we plot the Washington-DDO photometry along with a dividing line between dwarfs and 
giants (Paltoglou & Bell 1994 ). The boundary between giants and dwarfs is actually far 'fuzzier' than 
suggested by the solid line and complicated by the photometric transition from dwarfs to giants through 
subgiants. This soft boundary is readily apparent in Majewski et al. (2000) Fig. 14. Objects just above the 
heavy line are statistically more likely to be giants than objects just below the line. Reference stars DC-2 
and DC-5 have spectral types that are too early for this discriminant to work properly. The photometry is 
consistent with a giant or subgiant classification for the other reference stars. 



1 The Two Micron All Sky Survey is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/California Institute of Technology 

2 The Apache Point Observatory 3.5 m telescope is owned and operated by the Astrophysical Research 
Consortium. 
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3.3. Interstellar Extinction 

To determine interstellar extinction we first plot these stars on several color-color diagrams. A 
comparison of the relationships between spectral type and intrinsic color against measured colors provides 
an estimate of reddening. Figure || contains V-R vs V-K and V-I vs V-K color-color diagrams and reddening 
vectors for Ay = 1.0. Also plotted are mappings between spectral type and luminosity class V and III 



from Bessell fc Brett (1988) and Cox (2000, hereafter AQ00), again with reddening vectors and the loci of 
luminosity classes V and III stars. Figure along with the estimated spectral types, provides measures of 
the reddening for each reference star. 

Assuming an R = 3.1 galactic reddening law (Savage & Mathis 1977), we derive Ay values by 
comparing the measured colors (Tables |l| and ^) with intrinsic V-R, V-I, J-K, and V-K colors from Bessell 



fc Brett (1988)| and AQ00. Specifically, we estimate Ay from four different ratios, each derived from the 



Savage & Mathis (1977) reddening law: Ay/E(V-R) = 4.83; Ay/E(V-K) = 1.05; Ay/E(J-K) = 5.80; and 
Ay /E (V-I) = 2.26. These Ay are collected in Table [I| Colors and spectral types are inconsistent with a 
field-wide average <Ay > for the 6 Cep field. The spatial distribution of the average reddening star to star 
is shown in Figure [l]. A simple uniform extinction would predict a correlation between Ay and distance, 
with more distant objects having higher Ay. This correlation is absent in Table [| suggesting that either the 
extinction or the distances are in error. Alternatively, a patchy distribution of the ISM would destroy any 
correlation, a distinct possibility for this field at Galactic latitude, 1 = 0?5. As we shall see, the reddening 



for reference star DC-2 (discussed in Section 4.4), is of critical importance to an estimate of the reddening 
at the location of 5 Cep. 



3.4. Reference Frame Absolute Parallaxes 



We have prior knowledge that reference star DC-2 is thought to be physically associated with 5 
Cep. |Hoffleit fc Jaschck 1982| ( = BSC82) note common proper motion with S Cep and that DC-2 = 



ADS15987C (S Cep = ADS15987A), while de Zeeuw et al. (1999) include both 5 Cep and DC-2 in the 
newly discovered Cep OB6 association. Consequently, we first explored the astrometric properties of the 
four remaining reference stars, solving for parallax and proper motion of each in turn relative to a reference 
frame defined by the other three reference stars. We obtained no significant improvement in \ 2 by allowing 
any reference star (other than DC-2) to have a proper motion relative to the other three. In each case the 
spectrophotometric parallaxes discussed below entered the solution as observations. 
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We derive absolute parallaxes with My values from AQOO and the <Ay > obtained from the 
photometry. These are listed in Table ||, with three possible values for reference star DC-2, two depending 
on past spectral classifications. The last value for DC-2 is derived by constraining it to belong to the Cep 



OB6 association (see section [LJ below) . The weighted average absolute parallax for the reference frame is 
< T^abs >= 0.77 mas, including the highest weight parallax determination for DC-2, and 0.63 mas without 
DC-2. Statistically, DC-2 has very little weight in our reference frame. Nonetheless, it is astrometrically 



critical, as discussed in (Section 4.4), below. 



4. Absolute Parallax of 6 Cep 

4.1. The Astrometric Model 

With the positions measured by FGS 3 we determine the scale, rotation, and offset "plate constants" 
relative to an arbitrarily adopted constraint epoch (the so-called "master plate") for each observation set 
(the data acquired at each epoch) . The mJD of each observation set is listed in Table ||, along with the 
magnitude measured by the FGS (zero-point provided by Barnes et al. 1997), a phase (based on P = 



5.366316 days), and a B-V estimated by comparison with the UBV photometry of Barnes et al. (1997). 
The 5 Cep reference frame contains 5 stars. We employ the six parameter model discussed in McArthur et 
al. (2001) for those observations. For the 5 Cep field four of the reference stars are significantly redder than 
the science target and one is bluer. Hence, we apply the corrections for lateral color discussed in Benedict 
et al. (1999). 



As in all our previous astrometric analyses, we employ GaussFit ( Jcffcrys et al. 1987 ) to minimize x 2 . 
The solved equations of condition for 5 Cep are: 

x' = x + lcx(B - V) - AXFx (2) 

y' = y + lcy(B - V) - AXFy (3) 

£ = Ax' + By' + C + R x (x' 2 + y 12 ) - [i x At - P a7 r x (4) 

r) = -Bx' + Ay' + F + R y {x' 2 + y' 2 ) - fj, y At - P s n y (5) 

where x and y are the measured coordinates from HST; lex and Icy are the lateral color corrections from 
Benedict et al. 1999; and B — V are the B-V colors of each star, including the variable B-V of S Cep 
(Table |). Here AXFx and AXFy are the cross filter corrections in x and y, applied to the observations of 
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S Cep and reference star DC-2. 6 Cep has a full range of 0.2 < B-V < 0.6. For this analysis we linearly 
interpolate between the 1995 and 1998 cross filter calibrations (see Table 1, Benedict et al. 2002) as a 
function of 5 Cep color. A and B are scale and rotation plate constants, C and F are offsets; R x and R y 
are radial terms; fi x and fi y are proper motions; At is the epoch difference from the mean epoch; P a and 
P$ are parallax factors; and tt x and n y are the parallaxes in RA and Dec. We obtain the parallax factors 



from a JPL Earth orbit predictor (Standish 1990), upgraded to version DE405. Orientation to the sky is 
obtained from ground-based astrometry (USNO-A2.0 catalog, Monet 1998) with uncertainties in the field 
orientation ±0?05. 



4.2. Assessing Reference Frame Residuals 



The Optical Field Angle Distortion calibration (McArthur et al. 1997) reduces as-built HS T telescope 
and FGS 3 distortions with magnitude ~ 1" to below 2 mas over much of the FGS 3 field of regard. From 
histograms of the astrometric residuals (Figure ||) we conclude that we have obtained correction at the 
~ 1.5 mas level. The resulting reference frame 'catalog' in £ and rj standard coordinates (Table ^J) was 
determined with < er^ >= 0.3 and < >= 0.3 mas. 

Noting that the residual histograms have larger dispersions than we typically achieve, we plotted the 8 
Cep reference frame residuals against a number of spacecraft, instrumental, and astronomical parameters 
to determine if there might be unmodeled - but possibly correctable - systematic effects at the 1 mas 
level. The plots against residual included x and y position within the pickle; radial distance from the 
pickle center; reference star V magnitude and B-V color; and epoch of observation. Except for reference 
star DC-2 ( = HD 213307) discussed below, we saw no obvious trends, other than an expected increase in 
positional uncertainty with reference star magnitude. The largest residuals are associated with observations 
of reference stars DC-5 and DC-7 made during the two orbits with anomalous drift, discussed in Section 0. 



4.3. A New Companion for HD 213307? 

BSC82 notes a possible very short period companion (P<l d ) to reference star DC-2 = HD 213307. 
Such a companion would be undetectable by the FGS, either directly (changes in fringe structure) or 
indirectly (astrometric perturbation of the primary). Nonetheless, for DC-2 we found clear long-term and 
non-linear trends in the residuals with time. Because HST provides only relative proper motions, we do not 



- 9 - 



expect full agreement with the HIPPARCOS absolute proper motions. However, the agreement between 
HIPPARCOS and HST for 6 Cep was within the errors, while that for DC-2 was not (Table |). To assess 
the possibility that the residuals are caused by a perturbation due to a longer-period unseen companion we 
employed the model from our past binary star work (Benedict et al. 2001), 

£ = Ax' + By' + C + R x (x' 2 + y' 2 ) - fi x At - P a ir x - ORBIT x (6) 
i] = —Bx' + Ay' + F + R y (x' 2 + y' 2 ) - H y At - P s n y - ORBIT y (7) 

where ORBIT is a function of the traditional astrometric orbital elements. Due to the small number of 
epochs and observations we were unable to treat HD 213307 as we did Wolf 1062A (Benedict et al. 2001) 
and simultaneously solve for all the terms in equations 6 and 7. To obtain a solution we constrained all 
but the parallax, proper motion, and ORBIT parameters to values previously determined using equations 4 
and 5. The resulting orbital elements in Table |?j should be taken as highly preliminary. Figure ^| shows the 
residuals from the solution provided by equations 6 and 7 and the fitted orbit. Introducing these parameters 
to the modeling process reduced the residual histogram dispersions seen in Figure ^ from a x = 1.5 to 
(j x = 1.4 mas and from a y = 1.4 mas to a y = 1.0 mas. 

With M = 4Mq adopted as the mass of the B7-8 primary (AQ00), the perturbation orbit size, a = 2.0 
mas, and the period, P = 1.06 years yield a component B mass, Mb ~ 1.6Mq. If the companion is on the 
main sequence, this F0 V star would be approximately 7 mas (1.9 AU) distant from and ^3.8 mag fainter 
than the primary, consistent with the absence of a previous detection. Radial velocity variations of the 
primary would have an amplitude, Kl ~ 15 km s _1 , difficult to detect in a B7-8 III/IV star with vsini = 
135 km s- 1 (BSC82). 

4.4. The Absolute Parallax of 6 Cep and HD 213307 

In a quasi-Bayesian approach the lateral color and cross-filter calibration values were entered into the 
model as observations with associated errors. The reference star spectrophotometric absolute parallaxes also 
were input as observations with associated errors, not as hardwired quantities known to infinite precision. 
This approach allows us to incorporate any measurements relevant to our investigation. These include the 
HIPPARCOS parallaxes and proper motions of S Cep and HD 213307 (HIP 110988) with errors. 

Even though DC-2 is a binary with a poorly determined orbit, we find that we must include this 
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reference star in the solution. An ideal astrometric reference frame would surround 8 Cep. From Figure ^ 
it is clear that DC-2 is necessary to minimize extrapolation of the scale determined by the bulk of the 
offset reference frame. Without DC-2 the reference frame geometry is even less ideal. DC-2 also provides 
a constraint on the cross filter calibration, because both DC-2 and 8 Cep were observed with the neutral 
density filter. Following a suggestion from the referee, we constrained the difference in parallax between 8 



Cep and DC-2, using our prior knowledge of their proximity (see Section B.4). From dc Zccuw et al. 1999 



figure 25 we can estimate that the la dispersion in Galactic longitude for the OB association thought to 
contain both 8 Cep and DC-2 is 3°. One can therefore infer that the la dispersion in distance in this group 
is 3°/radian ~5%. Hence, the la dispersion in the parallax difference between two group members (e.g. 
DC-2 and 6 Cep) is 

Att = 5% X a/2 x 3.7 mas — 0.26 mas (8) 



where we have here adopted the mean parallax of Cep OB6, < n >= 3.7 mas, from dc Zeeuw et al. 1999 



This estimated difference allows us to assign a higher statistical weight to reference star DC-2 than if we 
had adopted either past spectral typing (lines 1 and 2 in Table |3|) or the HIPPARCOS value. The parallax 
difference between 8 Cep and DC-2 becomes an observation with associated error fed to our model, an 
observation used to estimate the parallax difference between the two stars, while solving for the parallax of 
<5 Cep. 

We obtain for 8 Cep an absolute parallax, 7r a fc s = 3.66 ± 0.15 mas, and for DC-2, ir a b s = 3.65 ± 0.15. 
Introducing the Cep OB6 parallax dispersion constraint and the HIPPARCOS parallax and proper motion 
measurements with their associated errors allows us to obtain a statistically significant result from a 
reference frame with very poor geometry. Our final 8 Cep parallax differs by ~ lamp and by ~ 3aHST from 
that measured by HIPPARCOS, ir a f, s = 3.32 ± 0.58 mas. Nordgren et al. (2002) have used long-baseline 
interfcrometry to measure the angular diameter of 8 Cep and other Cepheids. Their calibration of the 
Barnes-Evans (1976) relationship (surface brightness vs color index) yields a distance to 8 Cep of 272 ± 6 
pc and a parallax 7r a {, s = 3.68 ± 0.08 mas. We note that our 8 Cep absolute parallax and that of Nordgren 
et al. agree within each other's errors. 

In Figure || we compare 8 Cep astrometric parallax results from HST, HIPPARCOS, and Allegheny 
Observatory (AO). We plot the AO re-determination reported in Gatewood et al. (1998), not the original 
result (Gatewood et al. 1993). Also plotted is the Nordgren et al. (2002) result. Parallax and proper 
motion results from HST, HIPPARCOS, and AO are collected in Table ||. Parallax and proper motion 
results for reference star DC-2 are also presented in Table |8[ because it was measured by HIPPARCOS. 
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5. Discussion and Summary 

5.1. HST Parallax Accuracy 

Our parallax precision, an indication of our internal, random error, is often less than 0.4 mas. To 
assess our accuracy, or external error, we must compare our parallaxes with results from independent 
measurements. Following Gatewood, Kiewiet dc Jongc, fc Pcrsingcr (1998)| , we plot all parallaxes obtained 



by the HST Astrometry Science Team with FGS 3 against those obtained by HIPPARCOS. Data for these 
seven objects are collected in Table [| and shown in Figure 0. We have not included four Hyades stars whose 
parallaxes are considered preliminary (van Altena et al. 1997). The dashed line is a weighted regression 
that takes into account errors in both input data sets. 

The regression indicates a 2.5er scale difference between the HIPPARCOS and HST results. In addition, 
the fit x 2 suggests that either the HIPPARCOS or HST errors are overstated. Some, but not all of this 
discrepancy can be explained by pointing out that our log-log plot artificially spreads out these data which 
are effectively in two clumps, large and small parallaxes. Further exploring this issue, we have conducted 
a fully Bayesian analysis of the data in question, with an additional factor that represents the degree to 
which the presumed errors agree with the fit. We find that the input standard deviations of the HST and 
HIPPARCOS data points may have been overstated by a factor of ~ 1.5. It is not possible to decide which 
errors, HST or the HIPPARCOS, or some combination of the two, are overstated. However, HIPPARCOS 
errors have been subjected to many tests confirming their validity. Hence, it is likely that HST errors are 
overstated. 

Regarding the scale difference, we compare two hypotheses; the null hypothesis (a = 0, b = 1) and the 
complex hypothesis (a ^ 0, b ^ 1). We find that the Bayes factor against the null hypothesis (that the 
straight line in Figure ^ is correctly described by a = 0, b = 1) and in favor of a more complex hypothesis is 
about 30, depending on the priors we put on a and b. These 30:1 odds provide only modest support for the 
complex hypothesis, that there is a systematic scale deviation between the HST and HIPPARCOS results. 
We point out that the amount of data is very small, and that much of the scale difference (if real) depends 
upon the two largest parallaxes, which have a great deal of leverage. 

Measured proper motions provide another argument against the reality of this scale difference. Because 
it is desirable to reduce the impact of proper motion errors on an HST - HIPPARCOS comparison, we 
consider only two of the objects in Table ^| (Proxima Cen and Barnard's Star). They have proper motion 
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vector lengths exceeding 3800 mas yr -1 . HST yields proper motions relative to a local reference frame and 
HIPPARCOS produced absolute proper motions. The dissimilar approaches can result in unequal proper 
motions for the same object. However, comparing HST with HIPPARCOS, the average difference between 
these proper motion vectors is -0.01%, indicating a negligible scale difference. 



5.2. The Lutz-Kelker Bias 



When using a trigonometric parallax to estimate the absolute magnitude of a star, a correction should 
be made for the Lutz-Kelker (LK) bias ( Lutz fc Kelker 1973j ). Because of the galactic latitude and distance 



of 5 Cep, and the scale height of the stellar population of which it is a member, we use a uniform space 
density for calculating the LK bias. An LK algorithm modified by Hanson (1979) that includes the power 
law of the parent population is used. A correction of -0.015 ± 0.01 mag is derived for the LKH bias 
correction to the HST parallax of 5 Cep. The LKH bias is small because — ~ 4% is small. 



5.3. The Absolute Magnitude of 6 Cep 

To obtain the intrinsic absolute magnitude of S Cep we require the intensity-averaged apparent 
magnitude, the absolute parallax, and an estimate of interstellar extinction. To estimate the extinction we 
turn to reference star DC-2 (HD 213307). For this star we obtain an absolute parallax, 7r a f, s = 3.65 mas, 
between the two estimates depending on past spectral type and luminosity class determinations (Table |J) . 
From SIMBAD we collected B-V for all B8 III and B7 IV stars with V < 4 (presumably the nearest and 
least reddened members of their respective classes). The average for B8 III was <B-V> = -0.09 and for B7 
IV <B-V> = -0.13. We adopt an intermediate (B-V) = -0.11. From the measured color, B-V = -0.04, we 
obtain a color excess, E(B-V) = 0.07. Because a blue star can yield a different ratio of total to selective 
absorbtion, R, than a redder star, depending on the amount of color excess, we employ the formulation 
(Laney & Stobie 1993), 

R = 3.07 + 0.28 x (B - V)o + 0.04 x E{B - V) (9) 

and obtain for DC-2 < Ay > = 0.21 ± 0.03. The absolute magnitude of HD 213307 is then My = -1.10. 
Given that the distance to 8 Cep and DC-2 are the same within the errors, we adopt the same color excess, 
E(B-V) = 0.07 for <5 Cep, yielding (B-V) = 0.59 and R = 3.24. Therefore, < Ay > = 0.23 ± 0.03 for S 



- 13 - 



Cep. Our determination agrees with < Ay > = 0.25 ± 0.06 from Fernie et al. (1995)^. Because differential 
reddening effects are smaller at the I- band, we use a Savage & Mathis R = 3.1 reddening law to provide Aj 
= 0.14. 

We next require the intensity-averaged magnitude of 5 Cep. From Feast fc Catchpolc (1997)| we 
obtain an intensity weighted average <V> = 3.954. From the photometry of Moffett & Barnes (1985), 
we adopt a Johnson I-band intensity averaged magnitude, <I> = 2.993. With the HST absolute parallax 



from Section [L4j and < Ay > = 0.23 ± 0.03, we derive My = -3.47 ± 0.10, including the LKH bias and 
associated uncertainty. Similarly, Mj^qk) — ~ 4.14 ± 0.10, where we have transformed the Johnson to the 
Cousins-Kron photometric system as per Bessell (1979). The intensity-averaged photometry and absolute 
magnitudes are collected into Table [l(]. Using surface-brightness techniques, Gieren et al. (1993) obtain for 
5 Cep an absolute visual magnitude My = —3.63 ± 0.13 and, for log P = 0.73, My = —3.55 ± 0.11 from 
their Period-Luminosity (P-L) calibration for 100 Cepheids. Using a large number of relatively imprecise 
HIPPARCOS parallaxes, [Feast fc Catchpolc (1997)| and [Lanoix, Paturcl, fc Garnicr (1999] derive slightly 
different P-L calibrations for Cepheid variable stars. For log P = 0.73, Feast obtains My = -3.48 and 
Lanoix et al., My = -3.46. All determinations agree with our value within the errors. 



5.4. The Distance Modulus of the LMC 



Uncertainty in the distance to the Large Magellanic Cloud (LMC) contributes a substantial fraction of 



the uncertainty in the Hubble Constant (Mould et al. 200C). The HST Key Project on the Extragalactic 



Distance Scale ( Freedman et al. 2001 , Mould et al. 2000 ) and the Type la Supernovae Calibration Team 
( 3aha et al. 1999] ) have adopted the distance modulus value m-M = 18.5. Values from 18.1 to 18.8 are 
reported in the current literature, with those less than 18.5 supporting the short distance scale and those 
greater than 18.5, the long distance scale. Comprehensive reviews of the methods can be found in Carretta 



et al. (2000), Gibson (1999), and Cole (1998) 



Ideally, our absolute magnitude values for S Cep, My = — 3.47 ± 0.10 and M.j(ck) = — 4.14 ± 0.10, 
combined with apparent magnitudes for LMC Cepheids (corrected for reddening internal to the LMC) at 
log P = 0.73 would yield a unique LMC distance modulus. However, there are several complications that 
render such a distance modulus suspect. The first of these is cosmic dispersion in Cepheid properties due 



3 http://ddo. astro. utoronto.ca/cepheids. html 
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to the finite width of the instability strip. As shown in figure 5.13 of Binney & Merrifield (1998), lines of 
constant period are not horizontal in the HR diagram. This complication motivates the determination of 
Period-Luminosity-Color (P-L-C) relationships. Udalski et al. (1999) find a dispersion a = 0.07 magnitude 
about the P-L-C relation for LMC Cepheids. The scatter about any P-L relation decreases markedly 
as one utilizes progressively redder band passes that are less affected by internal reddening ( |Udalski et 



al. 199E and Lanvir 1999 ). A reddening- free magnitude, W, yields the P-L relationship with the least 
dispersion, Udalski et al. (1999) finding a = 0.08 mag. Assuming that the remaining scatter in a P-L-C 
or reddening-free P-L relationship is cosmic scatter, we adopt 0.075 mag as the cosmic dispersion in the 
intrinsic properties of Cepheid variables, including the defining member of the class, S Cep. 

Varying amounts of line blanketing due to intrinsic metallicity differences can also affect the apparent 
magnitude of a star. A comprehensive discussion of its effect on the Cepheid P-L relation is presented by 



Freedman et al. 2001. We adopt their correction, used for both V and I band and differential with respect 



to the LMC, of -0.2 ± 0.2 mag dex 1 . Specific to our calibration of the LMC P-L relation, for 6 Cep 



Andrievsky et al. 2002 find [O/H] = +0.01. The LMC has [O/H] - -0.4 (Kennicutt et al. 1998). Hence, if 



8 Cep had LMC metallicity, it would be 0.08 ± 0.08 mag brighter in V and I. 



We adopt the OGLE LMC V and I P-L relations (Udalski et al. 199£) because they are based on a very 
large number of Cepheids. Utilizing these to generate apparent, absorption-corrected Cepheid magnitudes 
at log P = 0.73, we obtain LMC V-band distance moduli m-M = 18.50 ± 0.13 or, corrected for metallicity, 
m-M = 18.58 ± 0.15. The corresponding I-band moduli are m-M = 18.53 ± 0.13 and m-M = 18.61 ± 
0.15, respectively. The errors include an RSS-ed 0.075 magnitude cosmic dispersion. We list these distance 



moduli in Table 11. Distance moduli corrected for metallicity are listed in the rightmost column as m-M = 
f(Z). At this stage in its maturity, this term introduces as much uncertainty as correction. 

All estimates agree within their respective errors with our recent determinations from HST astrometry 



of RR Lyr ( |Benedict et al. 2002[ ). There we obtained values of 18.53 ± 0.10 and 18.38 ± 0.10, where the 
range comes from differing <Vo(RR)> values for LMC RR Lyr variables found in the literature. We note 
that the errors associated with the Benedict et al. (2001) LMC distance moduli derived from RR Lyr 



neglected cosmic dispersion. According to Popowski & Gould (1998) this amounts to ~ 0.14 magnitudes, 
increasing our distance moduli errors to 0.18 magnitude. Table [l^ and Figure || summarize LMC distance 
modulus determinations based on HST astrometry and compare them with weighted averages of a number 
of results based on RR Lyr and Cepheid variables. A more comprehensive version of this figure, along 
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with the corresponding complete table and references, summarizing over 80 determinations based on 21 
independent methods, can be found on the web^. All LMC distance moduli based on HST astrometry 
are consistent with the LMC m-M = 18.50 ± 0.10 value adopted by the HST Distance Scale Key Project 
(|Frccdman et al. 2001]). 



5.5. Summary 

HST astrometry yields an absolute trigonometric parallax for S Cep, n a i, s — 3.66 ± 0.15 mas. This 
high-precision result requires an extremely small Lutz-Kelker bias correction, -0.015 ± 0.01 magnitude. 
To reduce our astrometric residuals to near-typical levels requires that we model reference star DC-2 as a 
binary and constrain it and S Cep to belong to the same stellar group, Cep OB6. Our astrometric results 
for DC-2 yield an extinction for that star of < Ay > = 0.21 ± 0.03. Correcting for color-dependent R 
differences, we find < Ay > = 0.23 ± 0.03 for 5 Cep. The dominant contributor to the error in the resulting 
absolute magnitude for S Cep, My = — 3.47± 0.10, remains the parallax. We find an LMC V-band distance 
modulus m-M = 18.50 ± 0.13, uncorrected for metallicity. This value is in agreement with our previous 
determinations with HST astrometry of RR Lyr and the value adopted by the HST Distance Scale Key 



Project (Freedman et al. 2001) 
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Table 1. Astrometric Reference Stars: Photometry 



ID 


V 


B-V 


V-R 


V-I 


V-K 


DC-2 a 


6.30±0.02 


-0.04±0.02 








DC-3 


13.47±0.02 


1.60±0.10 


1.10±0.04 


2.17±0.04 


4.22±0.06 


DC-4 


12.68±0.02 


1.83±0.10 


1.34±0.03 


2.56±0.03 


4.93±0.06 


DC-5 


13.68±0.02 


1.34±0.10 


0.96±0.04 


1.94±0.04 


3.69±0.06 


DC-7 


14.18±0.02 


1.60±0.10 


1.11±0.04 


2.22±0.04 


4.19±0.07 



Tutz fc Lute (1977) 



Table 2. Astrometric Reference Stars: Near-IR and Washington-DDO Photometry 



ID 


K 


J-H 


J-K 


M-T 2 


M-51 


DC-2 








-0.17±0.02 


0.02±0.02 


DC-3 


9.25±0.02 


0.76±0.02 


0.93±0.03 


2.22±0.02 


-0.07±0.02 


DC-4 


7.75±0.02 


0.93±0.02 


1.21±0.03 


2.58±0.01 


-0.1±0.01 


DC-5 


9.99±0.02 


0.71±0.02 


0.84±0.03 


1.97±0.02 


0.01±0.03 


DC-7 


9.99±0.02 


0.67±0.02 


0.88±0.03 


2.25±0.01 


-0.05±0.02 
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Table 3. Astrometric Reference Stars: Spectral Classifications and Spectrophotometric Parallaxes 



ID 



SpT 



V 



My 



-Kabs (mas) 



DC-2 
DC-2 
DC-2 
DC-3 
DC-7 
DC-4 
DC-5 



B7 IV a 
B8 IIP 
B7-8 III-IV d 
Kf III 
G8 III 
K3 III 
Gl III 



6.30 

6.30 

6.30 

13.47 

14.18 

12.68 

13.68 



-0.85±0.4 b 

-1.35±0.4 b 

-1.10±0.1 

0.6±0.4 

0.8±0.4 

0.3±0.4 

0.9±0.4 



0.28 
0.09 
0.21 
1.63 
2.05 
2.07 
2.13 



4.2±0.8 

3.2±0.6 

3.70±0.26 

0.6±0.1 

0.5±0.1 

0.9±0.2 

0.7±0.1 



Putz k Putz (1977) 



Wegner (2000) 



Savage, Massa, Meade, & Wesselius (1985) 



from membership in CepOB6 (de Zeeuw et al. 1999) 



Table 4. Astrometric Reference Stars: Ay from Spectrophotometry 



ID 


Ay(V-I) 


Ay(V-R) 


Ay(V-K) 


Ay(PK) 


< A v > a 


DC-3 


1.83 


1.43 


1.80 


1.45 


1.63±0.12 


DC-7 


2.35 


1.98 


2.13 


1.74 


2.05±0.15 


DC-4 


2.06 


1.90 


2.03 


2.32 


2.07±0.10 


DC-5 


2.37 


2.02 


1.96 


2.15 


2.13±0.10 


<Ay> b 


2.15 


1.83 


1.98 


1.91 


1.97 ± 0.13 



a average by star 
b average by color index 
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Table 5. S Cep Log of Observations 



Data Set 


mJD 


phase a 


V b 


B-V c 


1 


49908.60462 


0.870 


3.967 


0.63 


2 


49942.7928 


0.241 


3.910 


0.64 


3 


50087.07986 


0.128 


3.775 


0.55 


4 


50104.13513 


0.306 


3.989 


0.69 


■5 


50630.66157 


0.423 


4.103 


0.77 


6 


50668.40227 


0.456 


4.131 


0.79 


7 


50799.52478 


0.891 


3.862 


0.57 



Thase based onP = 5. d 366316, T = 43673.644 (mJD) (Barnes 



et al. (1997)) 



b Diffcrcntial FGS photometry. Zero point from Barnes et al. 
(1997). 

c Estimated from phase and Barnes et al. (1997) photometry. 
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Table 6. 6 Cep and Reference Stars: Astrometry 



ID 




71 a 






5 Cep 


45.2774±0.0003 


104.9195±0.0003 


0.0174±0.0002 


0.0050±0.0002 


DC-2 


37.4174±0.0003 


64.9353±0.0003 


0.0215±0.0003 


0.0037±0.0003 


DC-3 C 


0.0000±0.0003 


0.0000±0.0004 






DC-4 


131.7161±0.0002 


116.1802±0.0003 






DC-5 


88.9765±0.0003 


-27.4270±0.0003 






DC-6 


167.0683±0.0003 


76.5763±0.0003 







a £ and r\ are positions in arcseconds relative to DC-3 

\i x and n y are relative motions in arcsec yr _1 
C RA = 22 h 29 m 04.59 s Dec = +58°47' 40'.'7, J2000, epoch = mJD 50104.1363 



Tabic 7. HD 213307 - Elements of Perturbation Orbit 



Parameter 


Value 


a(mas) 


2.0± 0.2 


P(days) 


390 ± 9 


P (years) 


1.07 ± 0.1 


To 


2002.3 ± 0.1 


e 


0.35 ± 0.09 


i(°) 


36 ± 14 


n(°) 


100 ± 8 


O 


259 ± 7 
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Table 8. Final 8 Cep and HD 213307 (DC-2) Absolute Parallax and Relative Proper Motion Compared to 

Previous Results 



Parameter 5 Cep HD 213307 



HST study duration 2.44 y 

number of observation sets 7 

ref. stars < V > 12.06 

ref. stars < B - V > 1.3 

ref. stars < V > excl. DC-2 13.50 

ref. stars < B - V > excl. DC-2 1.6 



HST Absolute Parallax 3.66 ± 0.15 3.65 ± 0.15 mas 

HIPPARCOS Absolute Parallax 3.32 ± 0.58 3.43 ± 0.64 mas 

AO Absolute Parallax 2.8 ± 0.7 mas 



HST Relative Proper Motion 
in pos. angle 

HIPPARCOS Absolute Proper Motion 
in pos. angle 



17.4 ±0.7 21.8 ± 1.2 mas y" 1 

73°± 3° 80°± 5° 

16.9 ± 3.1 16.1 ± 2.7mas y" 1 

79°± 14° 74°± 12° 
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Table 9. HST and HIPPARCOS Absolute Parallaxes 



Object 


HST 


HIP 


HST Reference 


Prox Cen 


769.7 ± 0.3 mas 


772.33 ± 2.42 


Benedict et al. 


1999 


Barnard's Star 


545.5 ± 0.3 


549.3 ± 1.58 


Benedict et al. 


1999 


Feige 24 


14.6 ± 0.4 


13.44 ± 3.62 


Benedict et al. 


2000 


Wolf 1062 


98.0 ± 0.4 


98.56 ± 2.66 


Benedict et al. 


2001 


RR Lyr 


3.60 ± 0.20 


4.38 ± 0.59 


Benedict et al. 


2002 


8 Cep 


3.66 ± 0.15 


3.32 ± 0.58 


this paper 




HD 213307 


3.65 ± 0.15 


3.43 ± 0.64 


this paper 
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Table 10. S Cep Intensity- Weighted Magnitudes 



mag 


J a 


CK b 


DDO c 


< V > 


3.95 


3.95 


3.95 


<I> 


2.99 


3.18 




< V >-< I > 


0.96 


0.75 




Ay 


0.23 d 




0.25 


Aj 


0.14 d 






< V > 


3.73 


3.73 


3.70 


<Io> 


2.84 


3.06 




My 


-3.47±0.10 


-3.47±0.10 




Mi 


-4.34±0.10 


-4.14±0.10 





a Johnson VI photometry from Moffett & Barnes (1985) 
b Transformation per Bessel (1979). 



from http: / /ddo. astro. utoronto.ca/cephcids. html 



Estimated from Section 5.3 and Savage & Mathis (1977) reddening curve. 



Table 11. LMC Distance Moduli from the OGLE^epheid Period-Luminosity Relation 



Wavelength 



m-M b m-M = f(Z) 



b.c 



V 
I 



-2.760±0.031 17.04±0.021 15.03±0.03 18.50±0.13 18.58±0.15 
-2.962±0.021 16.56±0.014 14.40±0.02 18.53±0.13 18.61±0.15 



a Udalski ct al. (1999), m = a- log P + b 

b Errors are RSS of m error from PL, absolute magnitude (M) error, and 0.075 magnitude 
cosmic dispersion. 



Freedman et al. 2001 metallicity correction, A(m-M) = -0.2±0.2 mag dex 1 
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Tabic 12. V-band LMC Distance Moduli 



No. Object Source m-M 

1 S Cep this paper, m-M = f(Z) 18.58 ± 0.15 a 

2 RR Lyrae Benedict et al. 2002 18.53 ± 0.18 b 

3 all Cepheids Benedict et al. 2002, table 10 18.53 ± 0.07 c 

4 S Cep this paper, m-M^f(Z) 18.50 ± 0.13 d 

5 all RR Lyr Benedict et al. 2002, table 10 18.45 ± 0.08 e 

6 RR Lyrae Benedict et al. 2002 18.38 ± 0.18 f 



''LMC Cepheid <V> from Udalski et al. (1999); |Freedman et al. 2001 metallicity correction 



b Caretta et al. (2000) LMC <Vrr(0) > = 19.14 (includes +0.03, < [Fe/H] > correction). 
c Weighted average. Error is standard deviation of the mean from four independent techniques. 
d LMC Cepheid <V> from Udalski et al. (1999); no metallicity correction. 
c Weighted average. Error is standard deviation of the mean from five independent techniques. 
f Udalski et al. (1999) LMC <Vrr(0) > = 18.94 (+0.05, < [FeH] > correction) = 18.99 
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Fig. 1. — 5 Cep and astrometric reference stars. Symbol size is indicative of V mag nitude (Table p. The 



numbers within each identification box are the per-star <Ay > from Table ^, Section and E(B-V) per 
100 pc from Table |. 
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Fig. 2.— M-DD051 (M-51) vs M-T 2 color-color diagram for reference stars DC-3 through DC-7. The solid 
line is the division between luminosity class V and luminosity class III stars. Giants are above the line, 
dwarfs below. The reddening vector is for Ay = 1.0. The numbers are the reference star ID's plotted at the 
observed values. The circles are de-reddened values, based on the per-star <Ay > from Table 0. 
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Fig. 3. — V-R vs V-K and V-I vs V-K color-color diagrams for reference stars DC-3 through DC-7. The 
dashed line is the locus of dwarf (luminosity class V) stars of various spectral types; the dot-dashed line is 
for giants (luminosity class III). The reddening vector is for an Ay = 1. 
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Fig. 4. — Histograms of x and y residuals obtained from modeling <5 Cep and the astrometric reference stars 
with equations 4 and 5. Distributions are fit with Gaussians whose cr's are noted in the plots. 
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Fig. 5. — Absolute parallax determinations for S Cep. We compare astrometric results from HST, 
HIPPARCOS, and a recent determination from Allegheny Observatory (AO, Gatewood et al. 1998). 
Nordgren et al. (2002) have derived a parallax from a new calibration of the Barnes-Evans relation, denoted 
BE. 
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Fig. 6. — Reference star DC-2 (HD 213307) astrometric residuals from modeling with equations 6 and 7. 
Each observational epoch (denoted Data Set in Table |||) is coded with a unique symbol. The emphasized 
symbols are plotted on the derived orbit at the epoch of observation. The orbit is preliminary, and suggests 
a M = 1.6Mq companion. 
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Fig. 7. — Bottom: HST absolute parallax determinations compared with HIPPARCOS for the seven targets 
listed in Table ||. Top: The HIPPARCOS residuals to the error-weighted impartial regression line. The error 
bars on the residuals are HIPPARCOS la errors. 
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Fig. 8. — Recent determinations of the V-band distance modulus of the Large Magellanic Cloud (see 
Table [l2]). Values labeled 3 and 5 are weighted averages of a number of independent determinations using 
Cepheids (3) and RR Lyr variables (5). 



